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Exceptional points are singularities of eigenvalues and 
eigenvectors for complex vàlues of, say, an interaction pa- 
rameter. They occur universally and are square root branch 
point singularities of the eigenvalues in the vicinity of level 
repulsions. The intricate connection between the distribution 
of exceptional points and particular fluctuation properties of 
level spacing is discussed. The distribution of the exceptional 
points of the problem Ho + XHi is given for the situation of 
hard chaos. Theoretical predictions of local properties of ex- 
ceptional points have recently been confirmed experimentally. 
This relates to the specific topological structure of an excep- 
tional point as well as to the chiral properties of the wave 
functions associated with exceptional points. 



I. INTRODUCTION 

We address the qüestions: which are the common prop- 
erties of the quantum mechanical Hamilton operators 
that give rise to spectral properties that are ascribed to 
quantum chaos. If a matrix representation of a Hamilto- 
nian originating from a classically chaotic analogous case 
is given, we exploro the mathcmatical mechanism that 
produces the special statistical features of the spectrum 
for the particular parameter range where classical chaos 
is discerned. While classical chaos has no intrinsic sta- 
tistical property, a further puzzling question is: why can 
a plain statistical approach, with no physical input (such 
as GOE or GUE) (J, reproduce the statistical properties 
of the spectrum so successfully? 

We believe that the common root to the answer of these 
qüestions lies in what is called the exceptional points ||] 
(EP) of an operator. Most physical problems in quan- 
tum mechanics can be formulated by the Hamiltonian 
Hq + XHi where the parameter A can play the role of a 
perturbation parameter, or it may serve to effect a phase 
transition, or it may under variation steer the system 
from an ordered into a chaotic regime. The EPs of the 
full operator are the points A for which two eigenval- 
ues coalesce. Here we exclude genuine degenerades of 
the self-adjoint problem, in other words, the eigenvalues 
coincide for no real A. The EPs occur in the complex 
A-plane. Note that the operator is not self-adjoint for 
complex A- vàlues. 

The definition of EPs is general and applies also to 
operators in an infinite dimensional space, also when the 
spectrum of the operator has a continuum part. In the 
present work we restrict ourselves to finite ./V-dimensional 
matrices Hg and H\ as in this case the role of the EPs 
and the associated Riemann sheet structure is thoroughly 



understood (BflJ]. We do not believe that restriction to 
matrices has a major impact on our conclusions since 
virtually all the practical work even in connection with 
quantum chaos is donc in a finite dimensional matrix 
space. 

The physical significance of the EPs is due to their re- 
lation with avoided level crossing for real A-values. The 
spectrum Ek(\), k = 1, . . . , N has branch point singu- 
larities at the EPs, in fact, any pair of the N levels are 
generically connected by a square root branch point in 
the complex A-plane. If this happens near to the real 
A-axis, a level repulsion will occur for the two levels for 
real A-values. Globally, all the EPs determino the shape 
of the whole spectrum. There is a nice analogy to the 
more widely known connection between the singularities 
being poles of the scattering function and the shape of 
the cross section: in a similar way as the positions of the 
poles including their statistical properties determine the 
measurable cross section, the EPs determine the shape of 
the spectrum and in particular the oceurrences of avoided 
level crossings. The distribution of the EPs will therefore 
determine the fluctuation properties of level spacing. 

The positions of the EPs are fixed in the complex A- 
plane and are determined solely by Hq and H\ . For large 
matrices it is prohibitive to determine the positions of 
the EPs. However, it is possible to determine the dis- 
tribution reasonably well from the knowledge of the two 
operators. Generically, a high density of EPs is a suf- 
ficient prerequisite for the oceurrence of quantum chaos 
if they are randomly distributed according to a specific 
distribution function. 

In the following section we recapitulate the bàsics 
about EPs. Section three presents matrix models to ex- 
cmplify the distribution of EPs and level spacing fluc- 
tuations. Section four presents two recent experiments 
where (i) the local topological structure of an EP has 
been shown to be a physical reality and (ii) where the 
chiral property of the wave functions at the EP has been 
demonstrated. Section five presents a summary and Out- 
look. 



II. EXCEPTIONAL POINTS AND 
UNPERTURBED LINES 

Avoided level crossing is always associated with EPs 
HUH] wherever they occur for the levels £fe(A) of the 
Hamiltonian Ho + XH\ . We give an elementary example 
for illustration and briefly list the essential aspects with 
regard to EPs. 
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Consider a two dimensional matrix problem where Hq 
is diagonal with eigenvalues t\ and £2, while H\ is repre- 
sented in the form 



Hi = U ■ D ■ ZJ- 



(1) 



Here, the diagonal matrix D contains the eigenvalues u>\ 
and lü2 of the matrix H\ and U is the rotation 



U = 
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The eigenvalues of the problem Ho + XHi are 



(2) 



(3) 



where 
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Clearly, when 
lines 



the spectrum is given by the two 



Bg(A) = e k + \Lü k k = 1,2 



(5) 



which intersect at the point of degeneracy A = — (ei — 
£ 2)/(^i — W2). When the coupling between the two levels 
is turned on by switching on <f> the degeneracy is lifted 
and avoided level crossing occurs. Now the two levels 
coalesce in the complex A-plane where R vanishes which 
happens at the complex conjugate points 



Xr 



ei - e 2 

CJl — UJ2 



exp(±2i</ 



(6) 



At these points, the two levels Ek(X) are connected by a 
square root branch point, in fact the two levels are the 
vàlues of one analytic function on two diffcrcnt Ricmann 
sheets. 

These considerations carry over to an iV-dimensional 
problem Q. The diagonal matrix Hq contains the el- 
ements e k and D the elements uj k , k = 1,...,N; the 
matrix U is now an A^-dimensional rotation which can 
be parametrized by N(N — l)/2 angles. (In the quoted 
paper a parametrization was chosen so that U is unity 
when all angles are zero.) The EPs are determined by 
the simultaneous solution of the equations 



det{E -H - XHi, 







— det(E-H -XH 1 ) = 0. 
dE 



(7) 



There are generically N(N — 1) solutions which occur in 
complex conjugate pairs in the A-plane. At those points 
the N levels E k (X) are connected in pairs by square root 
branch points when they are analytically continued into 



the complex A-plane. Since the positions of the singular- 
ities determine the shape of the spectrum, and in partic- 
ular the fluctuation properties, a closer analysis is indi- 
cated. As is exemplified in the next section, the crucial 
condition for the occurrence of level statistics ascribed to 
quantum chaos is a high density of EPs in the complex 
plane within a small window of real A- vàlues. 

To get an idea about the density of EPs we introduce 
the concept of unperturbed lines. Clearly, when U is 
the unit matrix (all angles are zero), the spectrum of 
Hq + XHi = Ho + AD is given by the lines e k + Xuju with 
k = 1, . . . , N. The N(N - l)/2 intersection points of 
the N lines depend on the relative order of the numbers 
6fe and Wfe. If both sequences are in ascending order, 
all intersections occur at negative A-values; conversely, 
if one sequence is ascending and the other descending 
all intersections occur at positive A-values. In general, 
the order which is appropriate for the actual problem, is 
expected to lie between the two extremes. To find out 
the appropriate order we are guided by the asymptotic 
behavior of the levels Ek(X) of the full problem. For large 
vàlues of A the leading terms are given by 



E k (X) = Xu k + a k + 



(8) 



where the dots stand for first and higher order terms in 
1/A. Neglecting these terms, Eq.|| yields just the unper- 
turbed lines with the appropriate association of slopes Lu k 
and intercepts a k ■ From perturbation theory we find the 
latter to be the diagonal elements of the rotated Ho, viz. 

a k = (U- 1 ■ Ho ■ U) kM . (9) 



III. DISTRIBUTION FUNCTION 

We begin with the distribution function of the intersec- 
tion points of the unperturbed lines as this can be derived 
rigorously. Assuming a uniform random distribution for 
the u> k and the e k (the eigenvalues of H\ and Ho , respec- 
tively), the intersection points Xi ik — — (ej — efc)/(wj — u>k) 
are distributed according to @ 



P(x) 



const 



(10) 



The fact that P(x) cannot be normalized is due the fact 
that we assumed a uniform distribution of the uj k over 
the whole range of the real numbers; for a finite sample 
there will therefore be a deviation of P(x) around x = 
from the form given in Eq.(|Ï0|) as no \uj k \ is larger than 
some finite large number and henec an arbitrary small 
value of I/u is unlikely to occur. 

The intersection points are the points of degenerades 
of H + XH\ as long as U is equal to the identity matrix. 
We now gradually switch on the interaction between the 
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levels by switching on the angles randomly of the random 
orthogonal matrix U |tJ but keeping initially the interval 
from which the angles are chosen very small. In this way 
the degenerades become level repulsions as the EPs start 
moving out into the complex A-plane, a complex conju- 
gate pair from each degeneracy. We conjecture that their 
distribution function remains unchanged when they move 
out into the plane. For small vàlues of the mixing angles 
in U they will of course be concentrated around the real 
axis. The distribution function is now a function of two 
variables, say the real and imaginary part of a complex 
number. We parametrize this number by rexp(ia). For 
a fixed angle a the EPs are distributed according to 



P(r) 



const 



(11) 



This has been confirmed numerically in numerous cases. 
If the mixing angles of U are very small, there is an obvi- 
ous dependence on a as the EPs cluster around the real 
axis; yet for fixed a the distribution law is always given 
by Eq.(^ï]). Moreover, for the nearest neighbor distribu- 
tion (NND) of the energy levels it turns out that a proper 
Wigner distribution is obtained only when the EPs have 
fanned out into the plane in such a way that their dis- 
tribution becomes independent of the angle a. In Fig.l 
it is illustrated how the EPs fan out into the plane when 
the mixing angles of U are turned on. 



To summarize: hard chaos, that is a Wigner distribu- 
tion for the NND of the energy levels, is associated with a 
distribution of EPs in the complex plane that is indepen- 
dent of the angle a of the complex number r exp(ia) and 
depends only on the distance r according to Eq.([Í"H). If 
an a-dependence prevails, the NND is not Wigneresque. 
In the limiting case (U = identity) of the example used 
above the NND is a Poisson distribution. 




FIG. 2. Exceptional points as they move into the complex 
A-plane when the angles in U are switched on. Here a point 
of degeneracy at A = 5 is the starting point. In the figure on 
the right bottom this information is completely lost and hard 
chaos prevails. 

We only mention that in special non-generic cases the 
EPs are arranged in a geometrically ordered pattern likc 
in the integrable Lipkin model However, a small 
perturbation leads to the genèric situation as discussed 
above H. 



IV. LOCAL PROPERTIES OF EXCEPTIONAL 
POINTS 




FIG. 1. Exceptional points as they move into the complex 
A-plane when the angles in U are switched on. The lower part 
illustrates hard chaos. 

The law of Eq.([ïï]) is verified numerically even if the 
EPs originate out of a degeneracy of high multiplicity. In 
Fig.2 the emergence of EPs out of a point of degeneracy 
is illustrated. Here, Eq. (|ÏT|) is obtained once they fill the 
plane in an isotropic way, that is when the distribution 
has become independent of a. 



The topological structure of the square root branch 
point associated with an EP has been shown to be a 
physical reality in a recent experiment As a partic- 
ular consequence it has been established experimentally 
that the phases of the wave functions that take part in 
the coalescence of the two energy levels show a phase 
behavior that is distinctly different from that of a usual 
degeneracy, i.e. from that at a diabòlic point [fÜ)| . We 
cmphasise that one major signature, contrasting an EP 
from a degeneracy, is the lack of two independent eigen- 
functions; there is in fact only one eigenfunction: the EPs 
are the points where the Jordan form of the operator does 
not give a diagonal matrix J2jl. This is sometimes over- 
looked in the literature ||ÏT| , |l2|| . They can also occur in 
the continuum as the coalescence of two resonances [[Ï3) . 
The crucial experiment j^] yielded three major results 
which have been predicted in |Ï4| : 

1. If a loop is performed in the A-plane around the 
EP, the eigenenergies E\ and Ei are interchanged. 

2. The wave functions \tpi) and are interchanged 
by the loop and, in addition, one of them changes 
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sign. In other words, a loop in the A-plane trans- 
forms the pair {^1,^2} úito {—^2, V'l}- Therefore 
the two possible directions of looping yield difFerent 
phase behavior. In fact, encircling the EP a second 
time in the same direction, we obtain {—ipi, —1P2} 
whilc the next loop yields {^2, — Vh} and only the 
fourth loop restores the original pair. It follows that 
by going in the opposite direction, one finds after 
the first loop what is obtained after three loops 
in the former direction. This finding confirms the 
fourth root character of singularity for the wave 
functions. 

3. The eigenvalues E\,E2 have been studied as func- 
tions of A for two paths that were not closed. One 
path was just above, the other one just below A c . 
The results were different. Calling the real part of 
an eigenvalue the resonance energy and its imagi- 
nary part the resonance width. On one of the paths, 
the widths cross while the resonance energies avoid 
each other. On the other path, the resonance ener- 
gies cross while the widths avoid each other. 

We conclude the discussion with a furthcr important 
local property of EPs. The single and unique (up to 
a global factor) wave function at the EP has always a 
specific chiral behavior JÏ5| . A recent experiment at the 
TU-Darmstadt JÏ6| has directly confirmed this feature. 
In similar context this has been discussed for acoustic 
waves in a medium |Ï7| ] and indirectly observed in òptics 
[ IS | . The latter has been explained in terms of EPs in 
[ 1 9 1 . The essential finding of p5| is the unique form of 
the wave function at the EP being always of the form 



IV'EP) = IV'l) ± «l^) 



(12) 



where the plus or minus sign refers to a specific EP. The 
are the two wave functions that coalesce at the EP. 
No other superposition is possible, irrespective of a par- 
ticular physical situation such as driving the dissipating 
system. For a time dependent problem this signals chiral 
behavior. If the two wave functions relate to different 
parities or to different linear polarizations, the superpo- 
sition is obviously chiral; in the latter case it is a cir- 
cularly polarized wave of specific orientation. We recali 
that here Hq and H\ are assumed to be hermitian. If the 
two operators are non-hermitian there is still a unique 
superposition but Eq.(na) has to be modified p0|. 



V. SUMMARY AND OUTLOOK 

Exceptional points are a fascinating subject of thco- 
retical physics. As they are the only singularities of the 
spcctrum for a matrix problem, they 'make' the spec- 
trum, so to speak. They are directly associated with level 
repulsion. As a consequence, their statistical properties 



relate directly to that of the spectrum itself. Integrable 
systems give rise to a geometrical ordering while chaotic 
systems relate to disordered arrangements, yet with a 
specific distribution function. 

The local behavior is fascinating on its own. The 
topological structure of a square root branch point is a 
physical reality with all its consequences for the wave 
functions. In addition, with its intrinsic chiral behav- 
ior it may even hold some promise to shed light on 
the ubiquitous left-right asymmetry of our macroscopic 
world. On a speculative note: the kinematic relation 
E = ± yfp 1 + m 2 bears all properties of an EP at \p\ = 
±im. It was this relationship that led Dirac to his famous 
equation including spin and the properties of chirality. 



H.-D. Gràf, H.L. Harney, A. Heine, 
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